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This paper measures end-to-end (e2e) throughput of power beacon-assisted 
multi-hop decode-and-forward (DF) relaying scheme adopting 
non-orthogonal multiple access (NOMA) and transmit antenna selection 
(TAS)/selection combining (SC). Particularly, TAS/SC and NOMA are 
adopted at each hop to relay different data of a source to multiple 
destinations. Moreover, the transmitters including source and relays have to 
harvest wireless energy from radio frequency (RF) signals from a power 
beacon. We also propose a simple and efficient power allocation method for 
the signals transmitted at each hop. For performance measurement and 
comparison, we provide closed-form formulas of the e2e throughput over 
Rayleigh fading channel. We then verify our derivations by computer 
simulations as well as compare the e2e throughput performance between our 


Throughput 
Wirelessly energy harvesting 


scheme and the corresponding one that does not use NOMA. 
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1. INTRODUCTION 

Multi-hop relaying (MHR) [1]-[6] is commonly adopted to improve the end-to-end (e2e) 
performance for wireless networks when a source cannot directly send its data to a far destination. Therefore, 
MHR often employs intermediate nodes (called relays) to help the source-destination communication. 
In [1], [2], [5] the relays employed amplify-and-forward (AF) mode to forward the received data to their next 
node. In [3]-[6], the relays performed decode-and-forward (DF) mode, where they decoded, re-encoded, and 
then sent the received data to the next hop. The DF relaying often outperforms the corresponding AF one due 
to the noise removal. 

The main disadvantages of MHR are high e2e delay time and low throughput since the transmission 
is realized via multiple orthogonal time slots [1]-[6]. Indeed, if the number of hop is K, the maximum e2e 
throughput is 1/K, i.e., one data per K time slots. Although published works [4], [6]-[8] proposed various 
diversity-based MHR approaches (i.e., path selection [4], relay selection [6], transmit antenna selection 
(TAS)/selection combining (SC) [7], [8]) for performance enhancement, the obtained e2e throughput is still 1/K. 

Non-orthogonal multiple access (NOMA) [9]-[12] can be used to obtain higher e2e throughput for 
MHR. In NOMA, the source can simultaneously send different data to the destinations which can adopt 
successive interference cancellation (SIC) to extract the desired data. In [13], NOMA is applied in MHR 
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under presence of an active eavesdropper. Since the source can send two data to the destination at the same 

time, the e2e throughput obtained in [13] is 2/K. Moreover, the transmitters in [13] must adjust their transmit 

power, follows the instantaneous channel state information (CSI) of the eavesdropping links. 

Recently, radio frequency (RF) energy harvesting (RF-EH) [14], [15] has been reported in a lot of 
literature. One of efficient RF-EH models is to deploy power beacon stations [16]-[18] for charging the wireless 
devices within a certain area. Xu et al. [19], [20] studied the RF-EH based MHR employing the RF power 
beacon stations in underlay cognitive radio networks. Tin et al. [21] proposed the power beacon-aided MHR 
scheme, where NOMA is employed to send multiple data from the source to the destination. Vo et al. [22] 
studied the RF-EH NOMA MRR to optimize e2e secrecy performance. 

This paper proposes the RF-EH based DF MHR using power beacon and NOMA. In the proposed 
scheme, the source, relay and destination nodes are multi-antenna wireless devices, and TAS/SC is adopted at 
each hop to relay the source data. Next, the main difference between our proposed scheme and the existing 
ones related to MHR is presented as: 

— Different with [13], our scheme considers the RF-EH MHR method, and the transmitters and receivers 
in this paper are equipped with multiple antennas. Moreover, published work [13] studied physical-layer 
security with presence of the active eavesdropper. 

— Unlike [19], [20], the NOMA-based transmission is used to improve the e2e throughput, while the 
authors of [19], [20] considered the underlay cognitive environment. 

— The main difference between this paper and our previous work [21] can be listed as: i) this paper studies 
the multi-input multi-output (MIMO) scenario that employs TAS/SC; ii) the source in this paper 
attempts to send its data to multiple destinations while only one destination is considered in [21]; and 
iii) this paper does not consider the hardware impairments as in [21]. 

— The main difference between this paper and the related work [22] can be summarized as: i) this paper 
concerns with the e2e throughput performance; ii) the transmitter and receiver nodes in this paper are 
equipped with multiple antennas, and TAS/SC is adopted at each hop; and iii) this paper does not 
consider the physical-layer security as in [22]. 

Next, we will summarize the main contribution of this paper as: 

— We propose the wirelessly EH MHR scheme, where the TAS/SC and NOMA techniques are employed 
to obtain better performance, in terms of the e2e throughput and reliability of the data transmission. 

— We propose a simple and efficient power allocation method for the signals transmitted at each hop. 

— We derive closed-form expressions of the total e2e throughput for the proposed scheme over Rayleigh 
fading channel. We then verify our derivations via Monte Carlo simulations, i.e., the presented results 
show that the simulation and theoretical results are in an excellent agreement. 

— The proposed scheme can obtain better performance as compared with the corresponding MHR that 
does not use NOMA. 

The rest of this paper is organized as follows. The system model of the proposed RF-EH MHR with 
TAS/SC and NOMA is presented in section 2. Section 3 derives exact expressions of the e2e throughput. 
The simulation and theoretical results are given in section 4. Finally, section 5 concludes this paper. 


2. SYSTEM MODEL 

In Figure 1, the source Ty applies the NOMA technique to send N data to N destinations denoted 
Tkn, Where n E [1, N]. The transmitter Tẹ harvests RF energy from the single-antenna power beacon (B) for 
transmitting the data, where k € [0,K — 1], where B uses different frequencies with Tẹ. Without loss of 
generality, assume that the Tẹ and T,,, nodes have L receive and transmit antennas. In addition, the 
source-destination transmission is split into K orthogonal time slots, and TAS/SC is performed at each hop. 


o TAS/SC $ 


T, |sTAssc: 


Figure 1. System model of the proposed RF-EH MHR scheme with TAS/SC and NOMA 
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Let us denote Y,(a)y@) as the channel gain between the a-th antenna of the node X and the b-th 
antenna of the node Y. Note that in case that X or Y has only a single antenna, the indices a and b are 
skipped. When the channel is Rayleigh fading, y,(a),@) is an exponential random variable, and cumulative 
distribution function (CDF) and probability density function (PDF) of yy(a),(») are: 


Fy ayo (z) = 1 — exp(—Axyz) froyo (z) = Axy exp(—AxyZ) (1) 


Where: 
Axy = (dxy)* [21], B is a path-loss exponent, and dyy is the X-Y distance. 

If the e2e delay time is 1 (time unit), the time allocated for each time slot is tT = 1/K. Moreover, 
at each time slot, the time of at is spent for the EH phase, and the remaining time of (1 — @)T is spent for the 
data transmission. Considering the k-th time slot; the energy harvested by T;,_, is given as: 


Erg = natPg Emar Yor = natPg Xk- (2) 


Where: 
Pg is transmit power of B, n is energy conversion efficiency (0 < n < 1), and Xf) = Wee Yp: 


Then, the average transmit power of T,_, in the data transmission phase can be expressed as: 


ET Suoi sum 
Prag = Gage = APARE O 
Where 
u= ET After the EH phase, Tķ—ı linearly combines N signals to obtain a superimposed signal. 


xÈ = NP a Wa (4) 


Where: 
Xn is the signal of Tk n, an are coefficients: yy dy = 1and 1 >a, >a, >...> ay > 0. 
Before sending x , the nodes T;,_, and Tę cooperate to perform the TAS/SC technique as in [7]. 


Yoto = max, (mas (vro) ©) 


Where: 
tk-1 and rx are the selected transmit and receive antennas at T;,_, and Tk, respectively. 

After receiving xe from Tk-1, T; performs SIC to decode all the signals. Similar to [21], from (3), 
the instantaneous signal-to-noise ratio (SNR) obtained for decoding x,, can be written as: 


sum 
AnUAXk-1Y (tg) (rp) 
Tk-1 Tk 


Pr xn (n < N); andWx xy = An HAXRZIY (tp-1)..(rk) (6) 
Tk-1 Tk 


= N 5 
(Zinta Ti)HAXRIY (tp-1) (rg) tt 
Tr-1 Tk 


Where: 
A = P /0°, with o° is variance of Gaussian noises at Tẹ (also at all of the receivers). 
Then, the instantaneous channel capacity of x, at the k-th time slot is calculated as: 


Crx, = (1 — @)t log:(1 + Prz) = log, (1 + prx) (7) 


K 


Next, Tę repeats the operation that Tķ-1 did. Now, considering the last time slot (K-th time slot); 
Tx 1 sends the superimposed signal x®, = ei Pre 1 Xn to all the destinations. Also, the destination 
Tx+1n performs SIC to obtain its desired data x„. In addition, because the transmit power allocated to the 
signal xy is lowest, Tx, performs TAS for Ty y. Similar to (5), the channel gain of the Ty_, — Tx, link is 
written as: 
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= max | max b 8 
EER = 7Ra (o j 1ER) (8) 


Where: 
tx—1 and Tg y are the selected transmit and receive antennas at Ty_, and Tg,y, respectively. 
For the channel gains of the Tg-1 — Tx,» links (n + N); because SC is performed at Tk n, we can write: 


V pler- p(k) = PeDe CE (9) 


Where: 

Tx n 1s the selected receive antenna of Tx n- 

Similar to (6)-(7), the instantaneous SNRs and channel capacity obtained at the destinations can be given, 
respectively as: 


Han AX Y(t Koa) (Kn) 
K-1 Kn 
= - : n < N),and = WayAXe A 10 
Pk ,xn DI )HAXESIY (thea) ("e) t ( ) Wx xy Hay K 1Y plea) pH) ( ) 
(1-a) 
Cran => l0g2(1 + Wx xn) (11) 


Due to the DF MHR approach, the e2e instantaneous channel capacity of x, can be expressed by: 
Cercarn = MIN (Crx) (12) 


Finally, the total e2e throughput of the proposed scheme can be formulated as in [21]: 


(1-a) 
K 


TPNoma = Calin P (Coe. > Ca) (13) 
Where: 
Cy, is a pre-determined target rate. 

For baseline comparison, this paper also considers the RF-EH aided MHR without using NOMA 
(named OMA). In orthogonal multiple access (OMA), the e2e throughput of this model can be formulated as: 


(1-a) 
K 


Can Pr(COM^ > Ca) (14) 


TPoma = e2e 


Where: 


OMA _ x sum 
Cac = re ae (o — a)t logs (: + WAX 4-040) 


pe poii 


With 


= max | max b 
Vp (thea) plr) a=1.,.., (o iy ron )) 


3. PERFORMANCE ANALYSIS 
At first, we calculate Pr (Cozen > Cu) in (13) with n < N. Using (7), (10)-(12), we have 


sum 
Han AXk—-1Y (tk-1) (rk) 
Tk-1 Tk 


Pr(C > Ca) = TK, Pr(c > Ca) = [KP > On |, 
r( e2exXn = w) Hk=1 r( kxn = n) Hk=1 5 PREND Gi )MAXETLY (thea) (re) t! is 
k-1 k 


(15) 
Where: 

Cth 
On = 2070 — 1. In (15), if dn — On (Liner Ai) < 0, Pr(Corex, Z Cn) = 0. 
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Ortherwise, we rewrite (15) as: 


K 
Pr (Cezex 2 Cn) = pi Pr CERET 2 Pan) 
k-1 k 


k=1 
+00 ae 
Tha fy a rr) (> A ) fxg (z)dz ia 
Where: 
P = Oth 
ue [(an-Oin(renes a;))ua| 
In (16), Prad (.) is CDF of V pea) (4)? and from (5) and (9), we obtain: 


B edeo? z) = |z i r 1) o| = =(1 — exp(— Ar, “uo =1 ie Dch exp(— VAP, t,x)» ifk < K 


L 
7 vrv a 
P -dpl (z) = T o) =(1- exp(—Arg_sTknX)) =1+ yon C? exp( VArg_sTkn%) 


v=1 
(17) 
Also in (16), fgn (z) is PDF of X"). Because Xf = Yh, Vg using [7, eq. (41)], we have: 
k-1 
L 
(apr ) 2 
fxs sum m (z) = = S 1 exp(—Apr,_,Z) (18) 


Substituting (17) and (18) into (16), using [23, eq. (3.471.9)], we obtain: 


E: 2 
Pr (Corex 2 Cn) = Test Eze (L- oi ~(vAr,_ 17, ABT,_ Rey K,(2y vAr,_ 17, ABTp_ Pon)| 


2c? | 
rer oD bee TET Ce OR T (2 Parnes dota Pan) (19) 


Where: 
K,(.) is L-th modified Bessel function of the second kind [23]. 
With the same derivation manner, as n = N, we also obtain: 


K-1[ 1? 
Pr(Cosxy Z = Ca) = LI NG Lyte = n (vàr, 17, ABT, Pin)? 7K, (2 [Pare 17, ABTp_ Pun) 
k=1 |v= 
2 
2 (1 (L— D Eola 1TxnABTx- Pasa K, (2 [vrr 1Tgn^BTg- ‘Pus (20) 
Where: 
p — _ 4th 
BN — Canna 


Finally, substituting (19)-(20) into (13), we obtain an exact formula of TPNoma.- 
For the OMA scheme, similarly, it is straightforward to obtain that: 


a- -2 
TPoma = =” Ca M- [OHH pa (vir, 1TkABTp- ou K, (2 vAr,_,7,ABT,- oa)] 
(21) 
Where: 
— On 
Pn = Tay 
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Next, we propose the power allocation so that all the conditions ay — Om Yeins1 a; > 0 are satisfied. 
Indeed, we set an = 9n X? n41 di, Where € is a constant and ë > 1. Then, with Y“_, a, = 1, we have: 


_ Oh 2 1 
On = Tron” (asandan = AEO MNTT (22) 


Finally, we consider TPnoma and TPoma at high transmit SNR, i.e., A > +00. With (22), it is 


A>+0 
straightforward that PEC ate: > Cu) x 1(Vn). Therefore, we can approximate TPyoya and TPoya as: 


A>+0 (1-a)NC A>+ (1-a)C 
TPNomMa * ———* TPoma S — (23) 


The (23) implies that TPyoma can be higher N times than TPoma.- 


4. SIMULATION RESULTS 

Section 4 performs Matlab-based Monte-Carlo simulations [24], [25] to validate the derived fomulas 
of TPyoma and TPoya. We place the source (Tọ) and the relays (Tẹ) at positions(k/K,0), where 
k = 0,1,...,K — 1, and all the destinations are located at (1,0). In addition, the B station is fixed at (0.5, 0.5). 
For the illustration purpose only, the path-loss exponent, the target rate (C) and the energy conversion 
efficiency (7) are fixed by B = 3, Ca = 0.5 and ņ = 0.25, respectively. In Figure 2 to Figure 7 (all are 
drawn by Matlab), the simulation results are denoted by Sim, and the theoretical results are denoted by 
Theory. We can observe that all the simulation results match very well with the theoretical ones. 

Figure 2 presents the e2e throughput of the proposed NOMA scheme and the OMA one as a 
function of the transmit SNR A in dB, with different values of N. As we can see, the e2e throughput of the 
presented schemes increases as A increases because the transmit power of all the nodes increases. However, 
as A is high enough, the e2e throughput converges to the maximum value, as proved in (23). Moreover, as A 
is low, the performance of the proposed NOMA scheme is worse than that of the OMA scheme. Also, at low 
A values, the performance of our scheme is better with lower value of N. 

Figure 3 shows the e2e throughput of the NOMA and OMA schemes as a function of A in dB, with 
different values of L. Similar to Figure 2, the proposed scheme only outperforms the OMA scheme at 
medium and high A regions. It is shown that the e2e throughtput of both the schemes increases when the 
number of antennas (L) increases. However, at high A values, the e2e throughput of NOMA and OMA 
converges to the approximate values which do not depend on A. It is worth noting that the e2e throughput of 
both the schemes converges more rapidly with higher value of L. 


0.6 y 0.3 a 
Æ  OMA-Sim 
L OMA-Theory | 
0.5 @ NOMA-Sim (N=2) 0.25 
E NOMA-Sim (N=3) 
Vv NOMA-Sim (N=4) 
04+ NOMA-Theory 0.24 
5 + 
a 5 
S = 
D 0.3} © 0.15 oo 2 o oo o 
© = 
T g 
Fl ey 
a = Æ OMA-Sim (L=2) 
0.27 0.1 E NOMA-Sim (L=2) 
@ OMA-Sim (L=4) 
Vv NOMA-Sim (L=4) 
o1b 0.05 + OMA-Theory 
NOMA-Theory 
0 i i | i | | | 0 i | | 
5 0 5 10 15 20 25 30 35 -10 5 0 5 10 20 25 
A (dB) A (dB) 


Figure 2. Throughtput as a function of A (dB) when Figure 3. Throughtput as a function of A (dB) when 
K =3,€=2,a=0.1,andL = 2 K =3,¢ = 3, « = 0.1,and N = 2 
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Figure 4 investigates the impact of the number of hops (K) on the e2e throughput of the NOMA and 
OMA schemes. As observed, the e2e throughput of all the schemes decreases as K increases. We note that 
when K increases, the transmission time of each time slot decreases, which also decreases the data rate at 
each hop. Moreover, we see that with high value of K, the e2e throughput of the NOMA scheme is almost 
same when N changes. Therefore, our scheme does not perform well as the number of hops is high. 

Figure 5 studies the impact of € on the throughput performance. From (22), we see that the value of 
¢ is high, then a,_, is more higher than a,,. Therefore, this is the reason why as N > 2, the e2e throughput of 
the proposed scheme decreases as € increases. As ¢ is very high, the e2e throughput of our scheme is almost 
same. It is also seen that with N = 2, the performance of the proposed NOMA scheme only changes slightly. 
More interestingly, the proposed scheme obtains the highest throughput as N = 3. 

Figure 6 and Figure 7 investigate the impact of œ on the e2e throughput of the NOMA and OMA 
schemes. In Figure 6, we again see that the NOMA scheme is better than the OMA scheme as A is high. 
In both Figure 6 and Figure 7, it is worth noting that there exist optimal values of œ, at which the e2e 
throughput of the NOMA and OMA schemes is highest. However, as æ is very high, the e2e throughput of all 
the presented schemes decreases because the transmission time of the data transmission phase decreases. 
In Figure 7, the OMA scheme obtains better performance when a is small. Also in Figure 7, the performance 
of the NOMA scheme is worsest with N = 4, and with N = 3, the e2e throughput of the NOMA scheme is 
highest for a wide range of æ. 


2.5 T r r r T r r - 0.45 
i 
0.4 F 
ot 
0.35 + 
Æ OMA-Sim d iS 
@ NOMA-Sim (N=2) as 
15 E NOMA-Sim (N=3) 2 03 Ty 
a V _NOMA-Sim (N=6) ts 
5 OMA-Theory > a a E 
5 NOMA-Theory ° | Æ  OMA-Sim 
° 0.25 @ NOMA-Sim (N=2) 
E E NOMA-Sim (N=3) 
V NOMA-Sim (N=4) 
0.2 f x% NOMA-Sim (N=5) 
OMA-Theory 
NOMA-Theory 
0.15 J 
#1 E E E E E 
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Figure 4. Throughtput as a function of K when Figure 5. Throughtput as a function of € when 


A = 20 (dB), = 2, æ = 0.2, and L = 3 A = 10 (dB), K = 3, æ = 0.15, and L = 3 
0.3 r r r r 0.4 
0.35 
0.25 F 
0.3 F 
0.2 OMA-Sim (A = 0 dB) 
@ NOMA-Sim (A = 0 dB) 0.25 


%4  OMA-Sim (A = 15 dB) 
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Figure 6. Throughtput as a function of æ when Figure 7. Throughtput as a function of a when 
K=4,N =3,7 =2,andL =2 K = 2, A = 7.5 (dB), x = 2, and L = 2 
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5. CONCLUSION 

This paper proposed the wirelessly energy harvesting based multi-hop relaying scheme adopting the 
TAS/SC and NOMA techniques which can enhance the end-to-end throughput N times, as compared with the 
corresponding OMA scheme. The results also presented that the proposed scheme does not perform well at 
low transmit SNR regime and high number of hops. Moreover, the e2e throughput of the proposed scheme 
can be enhanced significantly by optimally designing the system parameters such as the number of hops (K), 
fraction of time allocated for the energy harvesting phase (æ) and the power allocation coefficient (€). 
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